Aspheric Surfaces

There is support for tracing over aspheric surfaces by specifying the base sphere radius of curvature at the vertex R, the conic constant K, and any optional, additional surface functions (e.g. polynomials) p.

This can be done either by passing keyword arguments for the latter two parameters to the real raytrace methods or by using the aspheric Layout constructor:

# example for a single surface, including object space
R = [Inf, -100.0]
t = [0.0, 0.0]
n = [1.0, -1.0] # reflector
K = [1.0, -0.7] # elliptical profile, prolate spheroid surface
p = [zero, y -> y ^ 4 / R[2] ^ 4] # fourth order surface sagittal correction
surfaces = Layout{Aspheric}(; R, t, n, K, p)

2×4 Layout{Aspheric}:
   Inf   0.0   1.0   1.0
 -100.0  0.0  -1.0  -0.7

There's also a Layout{Aspheric}(M) constructor which supports specification of the full 4 column matrix directly instead of using keyword arguments if the surface profiles are pure conic sections with no additional polynomial correction.